The two-loop supersymmetric corrections to lepton anomalous 
magnetic and electric dipole moments 

Tai-Fu Feng a ' 6 ' c , Xue-Qian Li d , Lin Lin a , Jukka Maalampi b , He-Shan Song a 

a Department of Physics, Dalian University of Technology, Dalian 116024, China 

~\q- b Department of Physics, 4OOI4 University of Jyvdskyld, Finland 

Q \ c Center for Theoretical Physics, Seoul National University, Seoul, 151-742, Korea and 

}_! ' Department of Physics, Nankai University, Tianjin 300071, China 

Q-i- 

<£ (Dated: February 2, 2008) 

^ ; Abstract 






*& 



cd 



Using the effective Lagrangian method, we analyze the electroweak corrections to the anomalous 
dipole moments of lepton from some special two-loop topological diagrams which are composed 
of neutralino (chargino) - slepton (sneutrino) in the minimal supersymmetric extension of the 



standard model (MSSM). Considering the translational invariance of the inner loop momenta and 



the electromagnetic gauge invariance, we get all dimension 6 operators and derive their coefficients. 



o 

Q-i! After applying equations of motion to the external leptons, the anomalous dipole moments of lepton 

are obtained. The numerical results imply that there is a parameter space where the two-loop 



supersymmetric corrections to the muon anomalous dipole moments may be significant. 
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I. INTRODUCTION 

At both aspects of experiment and theory, the magnetic dipole moment of lepton 
as well as the electric dipole moment draw great attention of physicists because of their 
obvious importance. The anomalous dipole moments of muon not only can be used for 
testing loop effect in the standard model (SM), but also provide a potential window to 
detect new physics beyond the SM. The current experimental world average of the muon 
magnetic dipole moment is [l| 

a ex P = 11 659 203 ± 8 x 1(T 10 . (1) 

Contributions to the muon magnetic dipole moment are generally divided into three 
sectors: QED loops, hadronic contributions as well as electroweak corrections. With the 
hadronic contributions which are driven from the most recent e + e~ data, we can get the 
following SM predictions J2j, |3j,|4J 

a SM = n 65Q 18QQ ± g q x 10 -10 

a SM = 11 659 175.6 ± 7.5 x 1(T 10 

af 1 = 11 659 179.4 ± 9.3 x 10 -10 . (2) 

The deviations between the above theoretical predictions and the experimental data are 
all approximately within error range of ~ 2a. Although this ~ 2a deviation cannot be 
regarded as a strong evidence for new physics, along with the experimental measurement 
precision and theoretical prediction accuracy being constantly improved, this deviation may 
turn more significant in near future. 

In fact, the current experimental precision (8 x 10 -10 ) already puts very restrictive bounds 



on new physics scenarios. In the SM, t 
amount to 19.5 x 1(T 10 and -4.4 x 1(T 10 



re electroweak one- and two-loop contributions 
respectively. Comparing with the standard elec- 
troweak corrections, the supersymmetric corrections are generally suppressed by A 2 /A 2 , 
where A EW denotes the electroweak energy scale and A NP denotes the supersymmetric energy 
scale. However, there is a parameter space where the one-loop supersymmetric corrections 



are comparable to that from the SM |6(. Since the one-loop contribution can be large, the 
two- loop supersymmetric corrections are possibly quite important 7(. 

Utilizing the heavy mass expansion approximation (HME) together with the correspond- 
ing projection operator method, the two-loop standard electroweak corrections to muon 



anomalous magnetic dipole moment (MDM) have been evaluated [8j|. Within the frame- 
work of CP conservation, the authors of Ref. [2| present the supersymmetric corrections 
from some special two-loop diagrams where a close chargino (neutralino) or scalar fermion 
loop is inserted into those two-Higgs-doublet one-loop diagrams. Ref. [l0( discusses the 
contributions to muon MDM from the effective vertices H ± W Zf, y, /i (-f^o)77 which induced 
by the scalar quarks of the third generation. 

In this work, we calculate the corrections from some special two-loop diagrams which 
are composed of internal neutralino (chargino) and (scalar) lepton lines. Since the electric 
dipole moment (EDM) of muon is also of special interest in both theoretical and experimental 



aspects 11], we as well present the lepton EDM here by keeping all possible CP violation 
phases. All the diagrams which we are going to calculate, were not discussed in literature. 
Besides, we first express our results in the form which explicitly satisfies the Ward identity 
requested by the QED gauge theory. In order to rationally predict the muon EDM, we 
certainly need to take the current upper experimental bounds on electron and neutron EDMs 
as rigorous constraints into account. Nevertheless, if we invoke a cancellation mechanism 
a,„on g differ s_et ri c contributions Q, or M e those sf_ of the & ,t 
generation to be heavy enough [14j . the loop inducing lepton and neutron EDMs restrict 
the argument of the /i parameter to be < 7r/(5tan/9), but no constraints on other explicit 
CP violation phases are enforced. 

Here, we apply the effective Lagrangian method to get the anomalous dipole momen- 
tums of lepton in this work. In concrete calculation, we assume that all external leptons 
as well as photon are off-shell, then expand the amplitude of corresponding triangle dia- 
grams according to the external momenta of leptons and photon. Using loop momentum 
translational invariance, we write the sum of the triangle diagrams which correspond to the 
corresponding self-energy in the form which explicitly satisfies the Ward identity required by 
the QED gauge symmetry. Then we can get all dimension six operators together with their 
coefficients. After applying the equations of motion for external leptons, higher dimensional 
operators, such as dimension eight operators, also contribute to the muon MDM and EDM 
in principle. However, the contributions of dimension eight operators contain an additional 
suppression factor m 2 /A 2 compared to that of dimension six operators, where m M is the 
mass of muon. Setting A NP ~ lOOGeV, this suppression factor is about 10~ 6 . Under current 
experimental precision, it implies that the contributions of all higher dimension operators 



(D > 8) can be neglected safely. 

We adopt the naive dimensional regularization with the anti-commuting 7 5 scheme, where 
there is no distinction between the first 4 dimensions and the remaining D — 4 dimensions. 
Since the bare effective Lagrangian contains the ultra-violet divergence which is induced by 



divergent sub-diagrams, we give the renormalized results in the MS scheme and on-mass- 
shell scheme jl5| respectively. The two-loop theoretical prediction certainly relies on our 
concrete choice of regularization scheme and renormalization scheme, however, our numerical 
results show that there is only tiny difference between the theoretical predictions by different 
regularization and renormalization schemes. We will discuss this problem in our other work. 
Through repeating the supersymmetric one-loop results, we introduce the effective La- 
grangian method and our notations in next section. We will demonstrate how to obtain the 
supersymmetric two-loop corrections to the lepton MDMs and EDMs in Section ILTTl In the 
Section HVl we study the dependence of the lepton MDMs and EDMs on the supersymmetry 
parameters numerically. Conclusions are presented in the last Section. 

II. OUR NOTATIONS AND THE SUPERSYMMETRIC ONE-LOOP RESULTS 

The lepton MDMs and EDMs that we will calculate can actually be expressed as the 
operators 

jC = 6 a la^l F 

C^^-'-dJa^^lF^. (3) 

Here, I denotes the lepton fermion, F^ v is the electromagnetic field strength, m l is the lepton 
mass and e represents the electric charge respectively. Note that the lepton here is on-shell. 
In fact, it is convenient to get the corrections of the loop-diagrams to lepton MDMs and 
EDMs in terms of the effective Lagrangian method, if the masses of all the internal lines are 
much heavier than the external lepton mass. Assuming external leptons as well as photon all 
are off-shell, we expand the amplitude of corresponding triangle diagrams according to the 
external momenta of leptons and photon. Then we can get all higher dimension operators 
together with their coefficients. As discussed in the introduction, it is enough to retain only 





Ei 



FIG. 1: The one-loop self energy diagrams which lead to the lepton MDMs and EDMs in MSSM, 
the corresponding triangle diagrams are obtained by attaching a photon in all possible ways to the 
internal particles. 



those dimension 6 operators in later calculations : 
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with T> ^ = d^ + ieA^ and uj t = (1 =F 7s) /2. At one-loop level, there are two triangle 
diagrams which contribute to lepton MDMs and EDMs (FigO}. After expanding the ampli- 
tude according to external momenta, the triangle diagrams determine following dimension 
6 operators together with their coefficients as 
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Here, A/", V, W denote the mixing matrices of neutralinos and charginos respectively, 
/, J = 1, 2, 3 are the indices of generations. We also adopt the shortcut notations: 
c w = cos# w , s w = sin6> w , Cg = cos/3, where # w is the Weinberg angle, and tan/5 = v 2 /v l is 
the ratio between the vacuum expectation values of two Higgs doublets. As for the mixing 
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matrices of sleptons and sneutrinos, R~, R D are: 

( (M 2 ) (M 2 



LL ' K LR/ \R. =m l 5.. , (i, j = l, ••-, 6) 



i2t(J^)i2 p = m ^, («, J = 1, 2, 3). (7) 

Those 3x3 matrices are denned as 

(Ml),, = (M 2 l ) u + m 2 /( 5 7J + m 2 ( S 2 - ±) cob2/W„ , 
(M 2 J 7J = (M 2 ) 7J + m 2 / rj - m 2 , 2 cos2/W„ , 

(M 2 ) 7J = (M 2 ), 7 + ~m 2 cos 2(3 (I, J = 1, 2, 3) , (8) 

where M 2 , M 2 , v4 e are the bilinear and trilinear soft breaking parameters in the lepton 
sector separately, and /z H denotes the /^-parameter in the soft supersymmetry breaking terms. 
Applying the equations of motion for leptons in Eq.fjSJ), we can get the lepton MDMs as 
. v o s. e 2 r d 4 q 
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with x % = m 2 /A 2 , and A NP denotes the new physics scale. The definitions of the functions 
Pi,2{xa/), ipi,2,3{x,y) can be found in appendix O If we change our notations to that of 
Ref. |6|, one can notice that those expressions are completely the same as the corresponding 
equations given in Ref. p. In order to obtain the above expressions, we have used the 
following identities which originate from the loop momentum translational invariance: 

d D q q 2 D f d D q 1 

. (10) 
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In the CP conservation framework, the supersymmetric one-loop contribution is approxi 
mately given by 

100GeV\ 2 



Aa^~13x 1(T 10 (— J tan/3, (11) 



when all supersymmetric masses are assumed to be equal to A NP , and tan/3 ^> 1. 

Correspondingly, the one loop supersymmetric contributions to the lepton EDMs can also 
be written as 
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Certainly, supersymmetry inducing operators Of 3 also contribute to the lepton MDMs. 
As we have seen above, only the operators Of contribute to the EDMs of lepton at one loop 
level. However, we will find that operators Of also contribute to both lepton MDMs and 
EDMs at the two loop order. 

III. THE TWO-LOOP SUPERSYMMETRIC CORRECTIONS 

In this sector, we analyze the two-loop supersymmetric corrections to lepton anoma- 
lous dipole moments. The two-loop supersymmetric corrections to the coefficients of those 

8 



operators in Eq. (J3J) originate from the two-loop self-energy diagrams of leptons, which are 
depicted in Fig. 121 The corresponding dipole moment diagrams are obtained by attaching 
a photon to these diagrams in all possible ways. In these diagrams there is no new suppres- 
sion factor, except a factor arising from loop integration, and the divergence caused by the 



sub-diagrams can be subtracted in the MS or on-shell schemes safely. It turns out that for 
some regions of the parameter space the two-loop results are comparable with the one-loop 
contributions J£J UJJ|. The reason for this is that the dependence of the two-loop results on 
the relevant parameters differs from that of the one-loop results. Among those two-loop 
contributions which have been analyzed in the literature, the corrections to muon MDM 
from the effective vertices j^Hq, ^ZHq induced by the scalar quarks of the third generation 
can be very well approximated by the formulaes 

Aa'< 2L ~ -0.013 x 10-10 m ^H ^n /3 

t H 

Aa h ; 2L ~ -0.0032 x i -io gV*»*gg!ff sign^/J , (13) 

b H 

where m~ and m- are the masses of the lighter t and b, A t b denote the trilinear soft breaking 
parameters of the t and b quarks, respectively, and M H is the mass of the heavy CP-even 
Higgs bosons. As for the two-loop diagrams where a close chargino (neutralino) loop is 
inserted into those two-Higgs-doublet one-loop diagrams, the correspondingly contribution 
can be approximated as 

A v it ■,, -,„ m/tan/3\ /100 GeV\ 2 . , . 

if all supersymmetric masses are set equal, i.e. /i H = m 2 = M A = A NP , and the U(l) gaugino 
mass m 1 relates to the SU(2) gaugino mass m 2 by the GUT relation m 1 = 5rra 2 /(3s 2 c 2 ) with 
the CP conservation assumption. Here M . is the mass of CP-odd neutral Higgs. Although 

nn 

other contributions in Ref . [9|, HPJ cannot be neglected also, they cannot be approximated 
as the succinct formulaes above. 

Among the two-loop supersymmetric diagrams under investigation, the corrections to 
the coefficients of operators in Eq. (j^J) originate from three types of graphs: the lepton self- 
energy diagrams, where there are two neutralinos and two sleptons; a chargino, a neutralino 
and two sleptons or two charginos and two sleptons as virtual particles in the loop (FigJ2J). In 
our previous works [16], we analyzed the contributions to the rare decay b — > s^y and neutron 




(a) 






FIG. 2: The two-loop self energy diagrams which lead to the lepton MDMs and EDMs in MSSM, 
the corresponding "triangle" diagrams are obtained by attaching a photon in all possible ways to 
the internal particles. 

EDM from the same topological two-loop diagrams which are composed of gluino, chargino 
(neutralino) , and squarks. Certainly, Figs|2] does not include all diagrams with internal 
slepton/neutrnalino (chargino) which contribute to the anomalous dipole moments of muon. 
Beside those diagrams in FigEl there are two-loop diagrams where a neutralino (chargino) 
one-loop self-energy composed of lepton-slepton or a slepton one-loop self-energy composed 
of lepton-neutralino (chargino) is inserted into those one-loop diagrams in Figd However, 
those diagrams belong to different topological classes, and we will analyze the corrections 
from those two-loop diagrams in our future works. We will adopt below a terminology 
where, for example, the " neutralino- neutralino contribution" means the sum of those triangle 
diagrams (indeed two triangles bound together), which have two neutralinos and two sleptons 
with a photon being attached in all possible ways to the internal lines. Because the sum 
of the "triangle" diagrams corresponding to each "self-energy" obviously respects the Ward 
identity requested by QED gauge symmetry, we can calculate the contributions of all the 
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" self-energies" separately. 

Since the two-loop analysis is more subtle than the analysis at one-loop level, we show 
here in some detail how to evaluate all the processes, which contribute at two-loop level to 
the theoretical prediction of the lepton MDMs and EDMs. Taking the same steps, which we 
did in our earlier works [16j . we obtain the following expressions for the relevant effective 
Lagrangian from the " neutralino-neutralino" self energy diagram: 
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with X> = ((q 2 — gi) 2 — m 2 )(q 2 — m 2 )(q 2 — m 2 )(g| — m 2 )(gf — m 2 ). Those tedious 

XaXp l J xl Bj E i Xa 

expressions of the form factors \.Af a ' b , ) (p = 1, •••, 4) and (Af c ' d , ) (p = 5, 6) are 

v x^xy p v x° a x { y p 

listed in appendix [Bj in order to express the sum of those corresponding triangle diagram 
amplitudes which satisfy the Ward identity required by the QED gauge symmetry explicitly, 
here we use the identities given in appendbJXl In a similar way, we can rigorously verify the 
following equations with those identities: 



f d D q x d D q 2 \ \x%h 
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FIG. 3: The diagrams which cancel the ultra-violet divergence of the two-loop diagrams, (g) rep- 
resents the counter terms which originate from the corresponding one-loop diagrams. Here, the 
corresponding triangle diagrams are obtained by attaching a photon to the internal charged slepton 
or chargino. 
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In fact, this is the direct consequence of the CPT invariance in the fundamental Lagrangian. 
The coefficients of high dimensional operators in Eq. (JT5|l contain ultra-violet divergence 
that is caused by those divergent sub-diagrams. In order to obtain physical predictions 
on lepton MDMs and EDMs, it is necessary to adopt a concrete renormalization scheme 
removing the ultra-violet divergence. In literature, the ultra-violet divergence is removed in 
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either the on-mass-shell renormalization scheme (a, |9| , or the simpler MS renormalization 
scheme 1 10, 14, 161. As an over-subtracted renormalization scheme, on-shell scheme looks 



more physical than MS renormalization scheme. However, there are at least two external 
legs of the divergent sub-diagram being the internal lines of the whole two-loop diagrams. 



This signifies that an artistic on-shell scheme is not more superior to the simpler MS scheme 
in our case. Certainly, theoretical predictions on lepton anomalous dipole moments depend 
on the concrete renormalization scheme. Here, we present firstly the renormalized results 



which are obtained in the MS renormalization scheme. We put the relatively complicated 
results which are obtained by the on-mass-shell scheme, in the appendix. 

For the two- loop " neutralino-neutralino" diagrams, the bare effective Lagrangian contains 
the following ultra-violet divergence 
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with e = 2 — D/2, where D denotes the dimension of time-space. Generally, the renormal- 
ization scale A RE and the new physics scale A NP should be of the same order in quantity, 
but there does not exist a compelling reason to make them equal. Here we keep the ratio 
x R = A 2 /A 2 as a free parameter in the expressions. In Eq. (|17|). we only retain the 
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operators Of that contribute to the lepton anomalous dipole moments, and the ellipsis 
represents the convergent parts of those coefficients. Certainly, the ultra-violet divergence 
contained in the amplitude of counter diagrams (the first two diagrams in Fig. EJ) will exactly 
cancel these divergences: 
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Adding Eq. (fTHjl and Eq. (fTHjl . we get the two- loop " neutralino-neutralino" corrections 
to the lepton MDMs: 
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-(x T x l) ) 1/2 Q Nd (x -x. ,x -x. ,x )Ref(?7 / )^(?7 J ).J?7 J ) t (£ 7 ) f 

V ;J x / jV,4V ;J ! Bi > y 0' B _. > y J \ \ Iff ) 3 p \ Iff ) ip \ Iff ) a j Vsjv / ai 

together with the lepton EDMs 



H X a f-j x'i 



9T v°v° 



2(4vr)2( SwCw )4 Ai> 



l 3 -,// ^iv,i l-^.j) ^6. 5 x v o i S e.) a ' v o J 



I J ' -B- ' v u ' E- ' v u 

• X Q ^J Xg 



Imfe / ).,(77 J ).,(77 J ) t (£ 7 ) f ) -Im(V).„(£ J ).,(£ J ) t (?/ ) f 

\^N'3(3\ lff/ip\ Iff ) ajV^ff) aij \\ <N ' J/3 V^jv /i/3 V'JV 'aj V ' N ' ai 

L X a Xp L CJ % X a ^j Xp 



x Im((£ / ). fl (£ J ).„(£ J ) t (f 7 )* ) -Im(V).,(?7 J ).,(?7 J ) t (n 7 )* 

V VS Af'J/3 V^AT/i/3 \>>ff'aj ^N'aiJ \\ 'JV /j/3 V 'jvA/3V 'N'aj^ < tf ' ai 

-x 1/2 n N Jx -x. ,x ;x- ,x )Im((ni) ,(^),(^ r ) t (£ 7 r ) f 
-x^O^^'X-.x 'x. ,x n )Imf(r7 7 ) ,,(r7 J ).Jr7 J ) t (^ 7 ) t 

JV,3\ ;J ' B . ' Y ! B > „0 / VV In']0\ »jv/«/3V /jv / aj^N' ai 

a. a L -^-a j Q 

The form factors are expressed as 

' 1 x a a 3 x p 

-\ (2 + lnxJL(x x )+ ft (x x )] +ip 1 (x x ) 



24 v 

H-^ 2 (x x )-- 

x /3 3 Z 



.X- 



a 3 ^3_ 

"£*2,2 l-£ > -^b / "T" X E . 03 £?2,2 I 3 ' J X e . / 



(9 3 X- 



to, "i 



s j (9 3 x . 



*£' ^ 



— I— P I rv> • rv> rv> • rv> rv> \ 

'*- N,l K^iJ ) *°e. ) "^^0 ) ^E. ) "^^O / 
' *i X a a ] Xp 

^ ^N,3 \ X ,J ? X e- ' "^v ' X B ■ ' ^yO ) 
1 "•% Xa a 3 x p 



11 1 a 

- \(2 + \nx Ii )p 2 (x x ) + (p3 ( x x )+ x e li2 ( x x ) 



2 *id 2 x~ E1 - : ' v "xS'~«i' 



— I— r IT* " T* T* * T* T* I 

1 ^l X a ^j Xp 

^^na{ X ,jj X F -1 X ^o 'i X e .i X ^o ) 

1 E > X a °i x» 



11 \ a 



where 



1 a l X a ^J Xa 



Qm,n\ X ll X 2J 



x m ln n x, - x" 1 ln n x 

1 1 2 z 



,(19) 



(20) 



(21) 



(22) 
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and other functions pi^, (fi,2,3, F 



N,i 



4) are defined in appendix O In Eq. (f2*T]) . 



all terms in the brackets are correct only for the naive dimensional regularization scheme 



and MS renormalization scheme. 

As for the two- loop "neutralino-chargino" corrections, we have 



o v ± 



Aa 2L '** 
i 1 



2\Fix T m , 
v i 1 i 1 



MiKUOlMUvih 



2{Atx) 2 s a c 2 I m c„ 

\ I w w ^ w 

X ^M,1^ ' X Si> X ^±'l X m.1 X v ) 

Xa 3 *-0 

+4x / (x x ,X x J^Re((Ai)t i (C c J )i J (^)L(^) 

X ^ A/,3 I" -T MA W X i>, ' X v ± 'l X El X v° ' 

-2(x T x + ) 1/2 Ref(A J ) t (C J ) f (£ 7 ) f (V 1 ) ■„) 

X ^M,3 W^' X v ±' X w„1 X v o) 

{2x iI x Y ,2 m hI 



30 



E i' *p 



m w c 



M(K)UOiMt&^ 



x ^M,4lU; X-, x , ;i. , x ) 



E j' *p 



as well as 



Ad 2L >x* 
i 1 



2x ,m 

i 1 1 



2(4vr)Vc 2 A, 



m c a 

w j3 



I —^lm((\ j y (C J ) f (v 1 )* (V 1 



N'jP 



F -F 

M,3 MA 



X 



(0;x D .,X x ±;x^,X xO J 

+2( Xj x +X ) 1/2 im((\i)i (cy (CV (d 

(0;x c .,x^;x^,x x J 

) i / 2 imf(A j ) t (c J ) f (e 7 ) f fa 2 

8 Xa ^J X^j 

(x ) 1//2 m 7 

x '-^((Ai)UOL(^)L(ei) 



N 1 30 



I N 1 30 



+ 



V2m w c fi 



J/3 



where the couplings are defined as 



E 3' *£ 



(O ta = ((^)/,(V) Ql - -7^— (**) (*-!)*(«),*) 



v^m 



w~/3 



(AiL = (^)/,((AT) laSw -(Ar) 2Q c w ). 



(23) 



(24) 



(25) 
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With the naive dimensional regularization scheme and MS renormalization scheme, the 
form factors are written as 



1 x a 3 x (3 

-7^\{ 2 + lnX n)[pl( x »^ x J-PlOV^.)] +<Pl(Xp.,X J 

LA I v ' L x a X /3 J J x a 



"^l(z ,Z- )- 



*£' E i 



d 3 



;?; 



a <9 3 a; 



"£2,2 I 3 ' ± ; x p 4 ) x p 



d 3 



b j <9 3 a; - 



"£?2,2 \ X „0 1 X P, .) 



X V E i 



+ 



F Af.l ' F M2 \ X ,J1 X £;. ' X V Q ' X e- ' ^v" ' ' 
' J ( D ? X a ^j X^ 



^m,3\^i ^e-) • r ,±) :r K.> 2- n J 



'J X f3 



if 1 <9 2 1 

-|2(2 + lnx R )^ 3 (^,^±) + V3(a: e .,X x ±)--^— <? 2 , 2 ( x x ± > x ^ ) } 



~i~r M5 {v', x P , x ,;x.,x ) , 



x« V x^ 



^M,4l^' X i>, ■' X +' ^k- ' X v J 



i *£ 



- i |(2 + lnx R )p 2 (x % ,x x J + ^ 3 (x % ,a; x jH--a;-.^-^ 2 (x x „,a; % ) 



+ ^ A/,6 l^) ^ ' X + 'l X B-'> X y0 



i X l 



(26) 



The tedious expressions of functions F M i (i = 1, • • • , 6) are put in appendixO With the 



naive dimensional regularization scheme and MS renormalization scheme, the resulting the- 
oretical predictions on lepton anomalous dipole moments from two-loop " chargino-chargino" 
diagrams are similarly expressed as 



±~± 



Aa 2L, X ±x 
ll 



X 



(4vr)V lI 



2m 2 



m 2 c 2 

w ff 



H^ohMU&Mt 



2m 2 

i J 

m 2 c 2 

w /3 



^((OM)M)UOl)] n oA^^,x x± -,x Dj ,x ) 



+4(x ± x 



± , ±. 



,1/2 



2 2 

rm m 



m 4 c 4 

w /3 



i Re((^),,(^) ii3 (^)l i (^)L 



+R*((^(^) v (OL(eJ)L 






J X. 



m 2 m , 

l J x ± 



■= X0 Kefir] 1 ).J7] J ).J V j y (f J )t )tt r Ax ,;x. ,x ± ;x. ,x . ) 

w /3 



V^' 



,(27) 
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and 



±^± 



i 1 



"(4vr) 2 <A NP 



K x ,i 



\V 2 < 



r m 

ij_ 

m 2 c 2 

w g 



H^oUO.MtKt) 



m 



m 2 c 2 a 



+2(x ± x ± 



1 777, 

'^-Im( (ri I ).Jr] J )Jr] J y (?/)M 

m 4 c 4 \ V 'Ch0\ 'cA/SV 'C'aj\ 'C'aiJ 

w /3 



rm 
4m 4 c 



m 2 m ± 
^ X ' iMijll) Avi) AV J J Xt) ] )V r Sx »x. ,x ,;x- ,x ,) 

W # ^ 



Here, 



v / 2m w c 



^ ^C,l V 2 '; J > ^j ' X ± ' ^B ■ ' X ± ) 



w /3 ^ 



— \ (2 + lnx R )U(a; x ± )+Pi(x.,x ± )1 +^i(x a; ± 



+^ 2 (x. ,X J - - 



5 3 



X 



d 3 x 



~Qi,i \X ± , X- ) -\- x 



d" 



d 3 x 



"<?2,2 \ X ± > "''S ) 



-4— n It "^ x " T T I 



J xi 



^ ^C,3 I 3 "; J ) -^P, ' X ± 1 "^S, ' X ± ) 



1 1 1 ijZ 

-l(2 + \nx R )^(x xr x Pi ) + -^(x xr x Dt ) + - ¥ — g2 ,(x xV x Pi ) 



— I— n IT ' T T * T* T i 



1 1 9 2 



_|_F ( rp • rp qr* • np rp \ 



(29) 



where the definitions of the functions F Ci , (i = 1, ■ • • , 4) can be found in appendix P 



Thus, we obtain the MDMs and EDMs of leptons in the MS renormalization scheme. 
However, the on-shell renormalization scheme is also adopted frequently to remove the ultra- 
violet divergence which appears in the radiative electroweak corrections [la]. As an over- 
subtract scheme, the counter terms include some finite terms which originate from those 
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renormalization conditions in the on-shell scheme beside the ultra-violet divergence to can- 
cel the corresponding ultra-violet divergence in amplitude. In the concrete calculation per- 
formed here, we need the following counter terms to cancel the ultra-violet divergence in the 
one-loop corrections to the vertex E*Xa ^ 



8C _ 



B*v u l 1 



V2& 



— SlJ + l^ZDjI UapS; 



,, i -(SZ^)ij5jj5 a /3 



+-^(SZ x o)i 3a S I j8ij 



{R\) jJ [M l ps v +M 2 ^ 



V2c 



p 



em 



\Flm s c a 

v w w /3 - 



'8e 8m 8m^ <5s w 8cA 

— -r : \Ou6ijbaf3 



m 



m 



+-(SZ[ J )j I 5i j 5ai3 + -(5Z\ : )i j 5 I j5 al 3 + -(SZ x o)/3 a SijSi 



(«t)i(3+J)M/9 ^_ 



V2e 



\6ij5ij5ai3 + -(8Z^) JI 8 ij 8 a p + -(SZ[) ij 5 IJ 6 af3 

e c,„ / 2 2 E 



1 



+-{8Z* )p a 8 IJ 8,, 



2" v 
em 



ij 



vt. 



j(3+J)N?p 



lJ 



\f2m s c a 

v w w p - 



' 8e 8m 8m 

— + — 1 - + — - 



m 



i j 



m 



s... 



8c r _ 



I SjjSaS, 



ij^a/3 



+7;(SZ[ i )j I 8i j 8 af 3 + -(SZDijSjjSap + -(8Z* )Ba^ijSij 

2 2 E 2 



(4)iJ^i K -( 3 °) 



Here, 8e represents the renormalization correction to electrical charge, 8m w and 8m 3 stand 
the renormalization corrections to the W-boson and lepton masses respectively, <5c w , 5s w 
as well as 8c„ are the renormalization corrections to parameters c w , s w and c , and 
(8Z t ' )ji, (8Z-)ij, (8Z x o) a/3 separately denote the wave function renormalization constants 
of leptons, sleptons, and neutralinos. In the on-shell scheme, we can fix those renormal- 
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20J. They include UV 



ization parameters by the mass-shell renormalization conditions 
divergence which cancel the corresponding UV divergence in amplitudes and the finite con- 
tributions which are determined by the on-shell condition, to the resultant expression of 
the finite amplitudes. In a similar way, we can write the counter terms for the one-loop 
corrections to the vertex v*Xa I 1 



8C 



!>*X a ll 



f \8ij8ij8ap + -(8Zl) ij 8 IJ 8 ali + -(8Z x +)p a 8 IJ 8 i 



+ 2^Z^)jl8ij8 a l3 



(Rt)M 



em 



iff 



[O. 



i,j 



\/2m s c a 

v w w (3 



' 8e 8m 



+ 
e m 



ij 
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H \SijSij5af3 + -{SZ^jrfijSap + -(8Zl) ij 8 IJ 5 a p 



m „ K c J 1 

1 



+ o($ Z x-)p^u$ij 



2 



(RDjjUfn L+ , (31) 



with (5Z.)ij, (5Z x -) a/ 3 are the wave function renormalization constants of sneutrino and 
chargino respectively. In order to shorten the length of text, we put the expressions of 
the theoretical predictions on the MDMs and EDMs of leptons in terms of the on-shell 
renormalization scheme in the appendix. 

So far, we have obtained all the corrections from two-loop super symmetric diagrams 
shown in Fig|21 Beside the two loop diagrams discussed here, it is well known that the two- 
loop Bar-Zee type diagrams also lead to significant contributions to the fermion MDMs and 



EDMs in the supersymmetric theory 14(. The Bar-Zee diagram corrections to muon MDM 



are discussed in JlQj , the contributions of Bar-Zee diagrams to muon EDM are analyzed in 
|l7|. Beside those two- loop diagrams that have been analyzed in literature and the diagrams 
presented in this work, there are still large amounts of two-loop diagrams that have concrete 
contributions to muon MDM and EDM. The present status of two-loop calculations cannot 
be considered as a complete analysis on MDM and EDM of muon in the framework of 
super symmetry. In the following section, we will only consider the corrections from those 
two-loop diagrams in Fig. Elto the one-loop supersymmetric theoretical predictions on lepton 
MDM and EDM through numerical method with some assumptions on the parameter space 
of MSSM. 

IV. NUMERICAL RESULTS 

With the theoretical formulation derived in previous sections, we numerically analyze 
the dependence of the muon MDM and EDM on the supersymmetric parameters in this 
section. Especially, we will present the dependence of the muon MDM and EDM on some 
supersymmetric CP phases in some detail here. Within three standard error deviations, 
the present experimental data can tolerate new physics corrections to the muon MDM as 
— 10 x 1CT 10 < Aa M < 52 x 10~ 10 . Since the scalar leptons v , /i^ 2 appear as the internal 
intermediate particles in the two-loop diagrams which are investigated in this work, the 
corrections of these diagrams will be suppressed strongly when slepton masses are much 
higher than the electroweak scale. To investigate if those diagrams can result in concrete 
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corrections to the muon MDM and EDM, we choose a suitable supersymmetric parameter 
region where the masses of the second generation sleptons are lying in the range M. < 1 TeV. 
In this work, we neglect all other possible sources of flavor violation except those due to 
the CKM matrix, and try to avoid ambiguities of the unification conditions of the soft- 
breaking parameters at the grand unification scale in the mSUGRA scheme. The MSSM 
Lagrangian contains several sources of CP violating phases: the phases of the [i parameter 
in the superpotential and the corresponding bilinear coupling of the soft breaking terms, 
three phases of the gaugino mass terms, and the phases of the trilinear sfermion Yukawa 
couplings in the soft Lagrangian. As we are not considering the spontaneous CP violation 
in this work, the CP phase of the soft bilinear coupling vanishes due to the tree level 
neutral Higgs tadpole condition. Moreover, for the model we employ here, the mass of the 
lightest Higgs boson sets a strong constraint on the parameter space of the new physics. As 
indicated in the literature [18], the CP violation would cause changes to the neutral-Higgs- 
quark coupling, neutral Higgs-gauge-boson coupling and self-coupling of Higgs boson. The 
present experimental lower bound on the mass of the lightest Higgs bosons is relaxed to 60 
GeV. In our numerical analysis we will take this constraint for the parameter space into 
account. 

As a cross check, we have compared our one-loop supersymmetric prediction on muon 
MDM in CP conservation framework with that obtained with corresponding Fortran sub- 



routine on muon MDM in the code FeynHiggs 21 1, and find a perfect agreement. In the 
two-loop sector, we check our Fortran subroutine on two-loop vacuum integrals with the cor- 
responding programs in the package FeynHiggs, and also find that they agree with each other 
very well. For guaranteeing validity of the results, we also independently develop certain 
programs for those two-loop integrals to check our two-loop integrals in Fortran code. 

Without losing too much generality, we will fix the following values for the supersymmetric 
parameters: M, = M. = \A.\ = 500 GeV , \m, I = \mJ = 300 GeV. Taking u„ = 
300 GeV, we plot the MDM and EDM of muon versus the CP phase 6. = axg(A.) for 
tan/3 = 5 or tan/3 = 20 in Fig0J As tan/? = 5, one-loop supersymmetric correction to 
the MDM of muon (Dash-Dot line) reaches 4 x 10~ 10 . With our choice for the parameter 
space, the two-loop supersymmetric corrections can approximately be as large as 30% of 
the one-loop results. For tan/3 = 20, one-loop supersymmetric prediction on the MDM 
of muon (Sold line) is about 17.4 x 10~ 10 , whereas the relative correction from two-loop 
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super symmetric contribution to one- loop result is approximated as 3%. In other words, the 
corrections of these two-loop diagrams turn more and more insignificantly along with the 
increase of tan j3. Actually, the third and fourth terms of neutralino-neutralino contribution 
(Eq. [T^jl and the third term of neutralino-chargino contribution (Eq. |2*3*|) dominate the 
corrections to muon MDM from the two-loop diagrams in Fig. |21 and those terms depend 
on the parameter tan/5 very mildly Because the CP phase 9- affects the anomalous dipole 
moments of muon through the mixing matrix R. for sleptons of second generation , this 
leads to that the variation of the muon MDM versus 9. is very gentle. Taking 9. = ±7r/2, 
the theoretical prediction on muon EDM approximates as 1.8 x 10~ 24 (e • cm) for tan/3 = 5, 
and 2.5 x 10~ 24 (e ■ cm) for tan/5 = 20. The muon EDM of this order can be detected 
hopefully in near future experiments with experimental precision 10 -24 (e • cm) [11]. 

Now, we analyze the variation of supersymmetric corrections to the muon anomalous 
dipole moments with the CP phase 9 X = arg(rra 1 ). Taking /i H = 300 GeV, we plot the MDM 
and EDM of muon versus the CP phase 9 1 for tan/? = 5 or tan/3 = 20 in FigEl Here, we 
find that muon MDM depends on the CP phase 9 l very gently. For tan/3 = 5, the one- loop 
correction to the muon MDM is about 4 x 10~ 10 . With our choice for the supersymmetric 
parameters, muon MDM is approximated as 5.5 x 10~ 10 when we include the corrections 
from those two-loop diagrams in Fig. El As tan/3 = 5, the one-loop contribution to the muon 
EDM originates from the "neutralino-slepton" diagram, and two-loop contribution mainly 
originates from the "neutralino-neutralino" diagrams. Because the concrete dependence of 
one-loop result on the CP phase 9 1 differs from that of two-loop result on 9 1 drastically, it 
is easy to understand why the correction to muon EDM from two-loop diagrams becomes 
dominant. When tan/3 = 20, the one-loop correction to muon MDM is enhanced, this leads 
to that Aa^ can reach 17.4 x 10~ 10 , the correction from those two-loop diagrams to muon 
MDM turn insignificant now. However, the theoretical prediction on muon EDM exceeds 
the precision of future experiment already. 

Perhaps the most interesting subject to study is the variation of muon MDM and EDM 
versus the CP phase 9 2 = arg(m 2 ). Taking yU H = 300 GeV, we plot the MDM and EDM 
of muon versus the CP phase 9 2 = arg(m 2 ) for large tan/3 = 5 or tan/? = 20 in Fig® 
Generally, the two-loop correction to the muon MDM is 30% approximately for tan/3 = 5. 
In the largest CP violation (9 2 = ±7r/2) case, the EDM of muon is large enough and can 
be experimentally tested with the experimental precision in near future: 10 -24 (e ■ cm). For 
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FIG. 4: The supersymmetric corrections to the MDM and EDM of muon vary with the CP violating 
phase 6. = arg(-A-) when // H = 300 GeV and tan/3 = 5 or tan/3 = 20, where the dash-dot lines 
stand for the results of one loop with tan/3 = 5, the dash-dot-dot lines stand for the results of two 



loop in MS scheme with tan/3 = 5, the short dash lines stand for the results of two loop in mass 
shell scheme with tan/3 = 5, the solid lines stand for the results of one loop with tan/3 = 20, the 



dash lines stand for the results of two loop in MS scheme with tan /3 = 20, and the dot lines stand 
for the results of two loop in mass shell scheme with tan j3 = 20. 



tan/3 = 20, the one-loop supersymmetric corrections to the MDM and EDM are enhanced 
drastically. Especially for the muon EDM, it reaches 2 x 1CT 22 (e-cm), which can be detected 
easily in the future experiment. 

In the above analysis, we always suppose /x H > 0. It is well known that the sign of the 
one-loop contribution depends on the relative sign of /x H and m 2 . Assuming CP conservation 
with 1 — 9 2 — 9- — 0, we plot the muon MDM versus /x H as tan/3 = 5 or tan/3 = 20 in 
Fig. The plot shows that the correction from the two-loop diagrams can be neglected 
safely when |/x H | < 100 GeV. Except the parameter /x H , another parameter tan/3 plays an 
important role in our analysis. Assuming CP conservation and setting /x H = ±300 GeV, we 
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FIG. 5: The supersymmetric corrections to the MDM and EDM of muon vary with the CP violating 
phase 1 = arg(m 1 ) when ^ H = 300 GeV and tan/3 = 5 or tan/3 = 20, where the dash-dot lines 
stand for the results of one loop with tan/3 = 5, the dash-dot-dot lines stand for the results of two 



loop in MS scheme with tan/3 = 5, the short dash lines stand for the results of two loop in mass 
shell scheme with tan/3 = 5, the solid lines stand for the results of one loop with tan/3 = 20, the 



dash lines stand for the results of two loop in MS scheme with tan /3 = 20, and the dot lines stand 
for the results of two loop in mass shell scheme with tan j3 = 20. 



plot the theoretical predictions on muon MDM versus tan/3 in Fig. |H1 Since the dominant 
contribution from the two-loop diagrams depends on tan /3 weakly, the variation of two- loop 
correction is not very obvious with the increase of tan/3. In other words, the correction 
from those two- loop diagrams turns insignificant in large tan/3 case because the one-loop 
supersymmetric prediction on muon MDM is proportional to tan /3. 

All the numerical results are obtained under the assumption M. < 1 (TeV). Since 
the supersymmetric corrections to the muon MDM and EDM contain a suppression factor 
A^ w /Ag , the contributions from the considered diagrams Fig El become insignificant if 
the supersymmetric scale A SUSY 3> A EW . 
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FIG. 6: The supersymmetric corrections to the MDM and EDM of muon vary with the CP violating 
phase 6 2 = arg(m 2 ) when fj, H = 300 GeV and tan/3 = 5 or tan/3 = 20, where the dash-dot lines 
stand for the results of one loop with tan/3 = 5, the dash-dot-dot lines stand for the results of two 



loop in MS scheme with tan/3 = 5, the short dash lines stand for the results of two loop in mass 
shell scheme with tan/3 = 5, the solid lines stand for the results of one loop with tan/3 = 20, the 



dash lines stand for the results of two loop in MS scheme with tan /3 = 20, and the dot lines stand 
for the results of two loop in mass shell scheme with tan (3 = 20. 



V. CONCLUSIONS 

In this work, we analyze some two-loop supersymmetric corrections to the anomalous 
dipole moments of muon by the effective Lagrangian method. In our calculation, we keep all 
dimension 6 operators. We remove the ultra-violet divergence caused by the divergent sub- 



diagrams in the MS and on-shell renormalization schemes respectively. After applying the 
equation of motion for muon, we derive the muon MDM and EDM. Numerically, we analyze 
the dependence of muon anomalous dipole moments on CP violation phases. There is an 
experimentally allowed supersymmetric parameter space where those two-loop corrections 
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FIG. 7: The supersymmetric corrections to the MDM and EDM of muon vary with /i H for tan (3 = 5 
or tan (5 = 20 in CP conservation framework, where the dash-dot line stands for the results of one 



loop with tan/3 = 5, the dash-dot-dot line stands for the results of two loop in MS scheme with 
tan/3 = 5, the short dash line stands for the results of two loop in mass shell scheme with tan/3 = 5, 
the solid line stands for the results of one loop with tan /3 = 20, the dash line stands for the results 



of two loop in MS scheme with tan /3 = 20, and the dot line stands for the results of two loop in 
mass shell scheme with tan/3 = 20. 



on the muon MDM are significant and cannot be neglected, meanwhile the EDM of muon 
can be large enough to be experimentally detected with the experimental precision of near 
future. 
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FIG. 8: The supersymmetric corrections to the MDM and EDM of muon vary with tan/3 for 
/i H = ±300 GeV in CP conservation framework, Where the dash-dot line stands for the results of 



one loop with /j, H = —300 GeV, the dash-dot-dot line stands for the results of two loop in MS 
scheme with /i H = —300 GeV, the short dash line stands for the results of two loop in mass shell 
scheme with [i H = —300 GeV, the solid line stands for the results of one loop with /x H = 300 GeV, 



the dash line stands for the results of two loop in MS scheme with /i H = 300 GeV, and the dot 
line stands for the results of two loop in mass shell scheme with fi H = 300 GeV. 
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